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NOTE 


reference  to  an  equation  by  a single  number,  thus 
(1>)  1 s reference  to  an  equation  of  that  number  in  the 
same  section  of  the  dissertation.  reference  to  equations 
in  other  sections  of  the  dissertation  is  made  by  giving 
the  section  number  followed  by  the  number  of  the  equation; 
for  example,  (2-5)  refers  to  equation  5 in  section  2. 
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PL*NE  STRAIN  IN  Li-MINiiTED  ORTHOTROPIC  STRUCTURES 


0.  Introduction.  It  is  well  ttnown  In  the  theory 
of  plane  strain  in  orthotropic  materials  that  if  F satis- 
fies the  equation 


3?* 


+ 


2(7/ 


0 


(1) 


where  a,  p,  f,  and  V depend  upon  the  elastic  constants  of 
the  materiel,  then  the  stresses  are  given  by 


^ f . = a2F  ^ _ agF, 

- by2>  m 7xy  axay 


(2) 


The  solution  must  also  satisfy  the  boundary  conditions. 

For  laminated  material  it  1s  also  required  that  the  stresses 
fyy,  -fxy,  and  the  displacements,  u,  v,  be  continuous  at  the 
layers  of  separation.  It  is  shown  in  this  paper  that  poly- 
nomial solutions  of  (1)  fail  to  meet  this  condition  except 
for  three  rather  trivial  cases.  Approximate  methods  based 
upon  the  concept  of  function  space  and  others  upon  replacing 
the  differential  equation  by  finite  equations  are  then  de- 
veloped. Several  examples  are  worked  out. 
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1.  Plane  strain.  In  problems  of  plane  strain, 
with  no  tody  forces,  the  stress  components  fxx,  T^y,  and 
fxv,  must  satisfy  the  equations  of  equilibrium 


and 


<}7xx 

= o 

ax 

ay 

Zhi  - o 

b X 

ay 

(l) 

Tz  X ~ 

Tyz  = o . 

(2) 

The  infinitesimal  displacements  are  denoted  by  u 
and  v (with  w = 0) . The  strain  components  are  defined  by 


xx 


: 

3x‘ 


e 


yy 


_ 3v 

- V 


-xy 


= fr( 


bu  + Ov' 
by  dx 


(3) 


with  exz  = eyZ  = ezz  = 0.  The  strain  components  must  sat- 
isfy the  compatibility  equation- 


^ 2exx  + d2eyy  __  gc)2exy 
by2  ^x2  ^x^y 


(4) 


For  orthotro^ic  materia.1  where  the  planes  of  elas- 
tic symmetry  are  ta«cen  as  the  coordinate  planes,  Hooce' s 
Law  may  be  written 


- 3 - 


©xx  = ^xxAx  - ^x^yy/1'  y * ^zx^zz/1^ 


eyy  = -tfxy'TxxAx  ♦ ryy/Ey  - ^zyfzzAz 
ezz  = -Oxz'J'xx/^x  - ^yz'Tyy/Ey  ♦ TZZ/EZ 


exy  - 'Txy/2/<xy 


(5) 


where  Ex,  Ey,  and  Ez , are  the  Young's  moduli  in  the  direc- 
tions of  the  respective  axes;  is  Poisson' e ratio,  the 

ratio  of  the  contraction  parallel  to  the  J-axis  to  the  ex- 
tension parallel  to  the  i-axis  associated  with  a tension 
parallel  to  the  i-exls;  and  /<xy  is  the  modulus  of  rigidity 
associated  with  the  x and  y axes. 


Since  ezz  = 0,  we  have 


"fzz  - ^xzTxx^/kx  + ^yzTyyEz/Ly  . 


(6) 


Substituting  in  (5)  we  set 


eyy  “ -^Txx  + ^Tyy 

exy  = 


(7) 


where 


^ = (1  - <rxz<rzx)/Ex 


~ (<7xy  + crxz°"zy)/:  x ~ ( ‘Tyx  + ^z^zxJ/^y 


y = (1  _ «?yZ<rZy)/Ey 


* = l/(2/^xy) 


(8) 
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and 

Solving  (7) 

where 


°lj/£l  - <7Jl/sJ* 
for  Txx>  t^y , and 

^xx  = 03  exx  + ^eyy 
'fyy  - ^exx  + ^yy 
7"xy  = 2/*xyexy 

X = a/(«r  _ £2) 

= ^/(aK  _ ^2) 

60  = V/(ar  - ^2) . 


(9) 

rxy  gives: 


(10) 


(11) 


If  a stress  function  F can  be  found  3uch  that 


'Txx 


32F 

dy2’ 


fyy  - 


a2F 

ax2’ 


a 2f 

3x3y 


( 12 ) 


then  the  equilibrium  equations  (1)  are  satisfied.  Substi- 
tuting in  (7)  gives 


;xx 


= tx&U.  - ^ &Ll 


3v2 


3x2 


eyy 


+ r 


a2F 

3x2 


exy 


3xdy 


(13) 


Substituting  (13)  in  the  compatibility  equation 

( 4 ) , 


gives 


Let  ff  = (7"/oOiy  and  K=  (V  - /S)/(cnr)^,  then  equa- 
tion (14)  becomes 


a4g  + 2 + = 0. 


2Z2  a^f 


dx  37} 


The  solutions  of  (15)  are  of  the  form 

F = R{F1(z1)  + F2(z2)i 


where 


21 

z2 


with 


x + i£xy  = x + l[f<  + (K2  - 1)^]^7} 
x + i£2y  = x + lC*  — (K2  - 1)^]^jj 

^ = (r/eO*fK  * (K2  - 1)*J* 

&2  = (r/o oty<  - (K2  - D*]*. 


(15) 


(16) 


(17) 


(18) 


In  terms  of  F,  the  stress  end  strain  components 


are  given  by 
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and 


where 


(20) 

(21) 


Since 


$i2ei 

♦ S t - 2vi!2  =»V  * Zfitf  *r  - 2 A2 

= ot^2  - 2{v  - p)&\.2  + ^ = 0 

we  have 

Si  = 5i2(2v  " 

(22) 

and  alao 

S^fe2  = *7* 

(23) 

$±2  ♦ S22  = 2 (V  - /*)/<*• 


(24) 
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2.  Stresses,  etralna,  c.nd  displacements.  Let 

pl  = (pn  ♦ inn)^int2 
n=0 

*2  = r (Qn  + iNn)z2n+2  . (1) 

n=0 


Then 

^4  = ro  E‘ti*p*2i(i*p*i)(tltp  <■  m„p)ip(kipj5lp*V 

^Iz  = E rI(£i-p»2)(it'-pn)(alctp  ♦ m1Ctp)iP('c*p^PxltjP. 

dz22  *=3  P-0 

(2) 


Substituting  In  (1-19)  gives 


XX 


= Z U+p+2)( fc+p+i)(£ip+2pK+p  + <52p+2Qiti-p 

U-o  p=e 

* + l,r2P*2Nttp)lP,'2(i*pJxl£yp) 


n 

= z 

K-0 


Z K ( -1 ) ^pi’1  ( K+P+2  ) ( ICi-p+1 ) (^P)  ( <5iP+2P 
- p=0 
p even 

<m-n. 


P 


* i2ll\.p)jp  • t t-D*^ x*-p»2)(x»p»i)  c;p) 

p=i 
p odd 


• (V^p  * <faPT2NK.p)yp]^ 


(3) 
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I 

ic-0 


+ 


Z (-l)ap(c*-pi-2)(  Kt-p+l)  (KtPJ^lPpt+p 

. P=0 

p even 

<s2p-icn>>jp  ♦ z',(-i)i»,tiitv2)(^i)(,;p; 

p-i 
p odd 


( <5^pM 


tc  +-p 


<^2pNt+p)yp 


(A) 


Txy  = £- 


K— 0 


Z (-l)^p(  K+pt-2)  ( k-i-p+1)  (£gpj  (^!P>1M 

Lp=o 

p even 


tc+p 


♦ ^2pnMn»p)Jp  * l:"'Vx>*p‘*U»p+2>U»P»i)(lt*p) 


P=1 
p odd 


• W1****  * <52?+1^,P)yp 


x*. 


(5) 


Substituting  in  (1-20)  gives 


3XX 


= E [1l  (-l)ip+1(^p-*-2)(t+p*-l)(k^p)  (Sip6iP 

k=o  Lp^o  k ' 


s+p 


p even 


♦ 52pe2V0)jf  - j;’  (-i)iP’i(i*p*2)U»pn)CtKP) 

P-1 
p odd 


* (Sipei“K*p  * ^2P£2Ntt+p)yP 


x*. 


(6) 
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m 

evy  = £ 
t=0 


2l'K  (-i)ip(ic»p>2)(^pn)(^p)(iipSifi-1.p 

. P-0 
p even 


♦ <$2pS2'Wjp  ♦ E V1>*p*i(lctp+2>(lc*ptl)  (K*p) 

P=1 
p odd 


(£lP^lUK+p  + <$2p^2NKi-p)yP 


(7) 


e = VT 
exy  xy* 


(8) 


Integrating  (6)  and  (7)  gives 


u = £ *£*  (-l)^P+1(kti-p+2)  OW1)  (<$ip£i*ic+l  + ^2P£2^K+p)yP 
tt-0LP=0 
3 even 

- 2I~Vl)fcp+*U*P+2)  (^l1)  ^lPclMk+P  +52P«2NKi-p)yPl 

p-=l  -J 

p add 


• xK+1  + f(y) 


(9) 


v = 


£ (-l)^p(ic*-p«‘2)(l£*p+1i)  (^ip5i2|,.+p  ♦ S2pS2«t*p)yptl 

K=0  LP-° 
p even 

♦ "£  Vl)^^*P+2)(K+P+1;  (S^SiM^p  ♦ ^2ps2N^p)yp+1 

P-1 

p ddd 


* xfc  + g(x) . 


(10) 
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tie  determine  f(y)  end  g(x)  ty  substituting  (9) 
and  (10)  in  the  third  of  the  equations  (1-3)  and  using 
(1-22)  and  (8).  This  gives 

g(x)  = 51  ( k+2)  + i‘2~1S2Nit)xls;+1  (11) 

K=0 

f(y)  = £ (-l)iP(p*2)(i1Ptl«1i(p  ♦ f2P,'1e2Np)yP,-l 

P—0 

p even 

+ Z?  (-l)^-^(p+2)U1P^1€1Pp  + ^2P+l€2Qp)yPtl-  U2) 
P=1 
p odd 
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3.  Laminated  construction.  7«'e  assume  that  the 
material  consists  of  three  layers  cemented  together.  The 
inner  layer,  or  core,  of  thictcness  2c  will  he  distinguished 
by  the  subscript  1,  and  the  two  outer  layers,  or  faces, 
which  are  assumed  to  be  similar  and  of  thickness  b - c 
by  the  subscript  2.  Thus  Exl  is  the  Young's  modulus  of 
the  core  in  the  x-direction. 

At  the  planes  of  separation,  y = ±c,  it  is  neces- 
sary that  ^yy»  TXy , u,  and  v be  continuous.  Equations  ex- 
pressing this  condition  are  obtained  from  equations  (2-4), 
(2-5) , (2-9),  and  (2-10)  by  adding  and  then  subtracting  the 
value  for  y = c from  that  for  y = -c  and  then  setting  the 
coefficients  of  each  power  of  x equal  to  zero.  The  results 
are  as  follows: 

4m- 1 

Z (-l)^p(  tt+p+2)(  f-p+1)  cP(S11pPt+p>  i + S21PQfc*p,l 
P-0 

p even 

-fl.2Pplc+p,2  ~ ^22PQic+p,  2 ) = 0 (!) 

Q| 

"z  ( -1)  ^p“^(  K+p+2)  ( fc+p+1)  (K;PJcp(Snp+1Pfe,p,i 
P=1 
p odd 


+ 52iP^1Qtc+p,i  - Si2P+lpic+-p,2  - 522P+1Qk>p,2)  =0  (2) 


12 


13  ( -1)  ^p(  K«-p*2)  ( t?!1)  cP(^llP<£llpc+p,  1 + S2ip€2i(iICi.p>i 

P~0 

p even 


- <5'i2p£12Ptti-p,2  “ ^22Pfe22^li+p,  2 ) ~ 0 ^ 


^m- k.  i 

£ (-l)-p(  K+pt2)  (t+p'’1)  cPC^uPSuPk+p,  1 + ^21PS2iQk^p,1 

p-0 

p even 


- <Si2P^i2Pic+p,2  ” ^22P^22^ke+p,  2 ) = 0 


(4) 


23  (-l)^P  ^(p+2)cp(<S11p  1€11Pp>1  + ^21P  l£2r*p,l 

P=1 
p odd 


- *12P+l€12*p,2  - *22P^l622£p,2)  = 0 


(5) 


23  (-l)^P  ^(iti-p+2)  (n+p+1)  (*  Pj  cp(5nPMttfPji  + 52ipNc+p,l 

p=l 
p odd 


- <T12PMKl.p>2  - <S22PN fc+p,  2 ) = 0 


(6) 


ZQ  (-1)^P(  K+p+2)(  k+p+1)  (£^pJcP(illpi'1M£+p>l 

p even 


* - <5l2Pn“^p,2  - S22PnNSfp>2)  =0  (7) 
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Z (Kt?i1)cp(S11p£l1Mt>p,l  + 52lp€21Ntt+p,l 

p odd 

- S12Pe12^tti-p,2  - *22Pc22i;lt+p,2)  = 0 (8) 

Z (-l)£pt^(  it+pi-2)  l)  cprl(SiiP5iiMK;+pj  i 

P=-l 
p odd 

T ^21P^21Nk;-»-p,  1 - £l2P£  12^kc+p,  2 - ^22Pj22Nlc^p,  2)  = 0 

(9) 

(-l)^P+^(P+2)cp(511PS11Mp(1  + ^21PS21Npf  1 

P=1 
p odd 

- <Si2PSl2up,2  “ ^22P522lNp,2)  = 0 (10) 

In  equations  (1)  - (4),  (6),  and  (7)  the  range  of  k is 
given  by  0 £ & $ m,  in  equation  (3)  by  -1  $ k $ m-1, 
and  in  equation  (9)  by  1 £ «t  £ m+1. 

If  m > 0,  then  putting  k = a in  (1),  (3),  and  (4), 
and  k;  = e-1  in  (2),  we  get 

pml  + Qml  ~ pra2  - ^m2  = 0 
^ll2pal  * ^212(^ral  “ $l22pm2  “ ^222^ra2  = 0 
*llpml  ♦ £21^ml  - 612pm2  - e22^m2  = 0 
Sll^ml  + ^21^ml  ~ ^12pm2  - 522^2  = 0 UD 
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In  order  for  these  to  have  a non-zero  solution,  It 
Is  necessary  that 


1 

-1 

-1 

2 

11 

^212 

i 

on 
t— ' 

ro 

ro 

OJ 

CM 

<*> 

11 

*21 

i 

<*\ 

ro 

*€22 

11 

S21 

“512 

-S22 

- ^212_^ll2^S222“  (Qt2”Cil  ^ ^2"  ^1  il  “ 0 

(12) 

Since 

h*  ' hi2  = -2(ri/<*i>*<Ki2  - U* 

the  condition  for  non-zero  roots  reads: 

At  least  one  of  the  following  relations  must  hold 


(/®2  '^l)2  = (a2  ” oti)(^2  - /i).  (13) 

For  m = 0 equation  (2)  does  not  apply  and  we  have 
the  equations 
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P01  + Q01  - p02  “ ^02  = 0 
€11P01  + ^1^01  “ C12f02  - €22<<02  = 0 
^ 11^01  + 5 21^01  “ $12*02  ” $22^02  = 0 (!*) 


Their  solution  is 


101  = (62I-£22)  (^22-^12)  - ( e22~€12)  (^21- s22)qog 

(621"€ll)  ($11-$12)  ~ (ell-e12)  (i21*Sll) 

0 = (£22-€ll)  (Sll-Si2)  - (ell-£j2)  ( S22-sll)n 

01  ( €21~£Il)  ^11*^12)  * ^eil-ei2)  ^2l'Sll)  02 

1 _ _ (£21-€ll)  (^11-^22)  - (£ll-g22)  (S21-Sll)n 

^e21-eil)^ll-512)  - (Cll-C12)  (S21-Sll)  02 


(15) 

rutting  k = m-1  in  (6)  and  (3),  «c  = m in  (7),  and 
a = in-t-1  in  (9)  gives,  for  m > 0, 

*llMml  * 52lNml  ' S12Mm2  - *22Nra2  = 0 
^ll^llSl  + ^21£21Nml  “ 512£12Mm2  * 622£22Nm2  = 0 

^ll1^  11mid1  + 52l1^21Nml  - ^12l€12Mm2  - ^22lc22Nm2  = 0- 

(16) 

These  have  the  solution 


« ,=_S22(£  22-£12)m  _ n ,-^22(e22-^12)M  u - & 2m 

ml  1 — n 6 — i «2,  ml  7 — 72 7i-rN«2*  km2“- r-J‘m2* 

®11'*21~  ll'  <y2l(£2l-£ll)  S\2 


(17) 
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Putting  k - ra-3  in  (6)  and  t = m-2  In  (7)  gives, 
for  m 7 2 

-2,1  * *21wm-2,l  ~ *12Mm-2,2  “ *22iNm-2,2) 

= (mi-2)(m+l)(S115u;n>1  + ^2l5Nm,l  " *123Mm,2  " 5225Nm,2) 

2^llMm-2, 1 + ^21^111-2,1  " <S12Mm-2,2  “ ^22Nm-2,2) 

= (m+2)(m+l)(i1^  Mm>1  *■  *215Nm,l  “ *125lsm,2  " 5225nh»,2)- 

(18) 

These  are  Incompatible  with  equations  (16).  Therefore  the 
only  solutions  to  be  considered  are  for  m = 0,  1,  2. 

For  a = 0 we  have  only  the  equations  (16)  which 
have  the  solution  (17). 

For  m = 1 we  have  equations  (16)  and  In  addition 
the  following  equation  obtained  from  (10): 

^11S11M11  * 6215  21^11  - Sl25l2M12  ~ *22S22N12  = 0.(19) 

In  order  for  these  to  have  a non-zero  solution,  we  must 
have 


£ n 

&21 

*12 

. *^22  ' 

= 0 . 

£nfin 

*21*21 

-^12*12 

“ *22*22 

(20) 

hllzn 

621l521 

■612l512 

- i22l522 

snsn 

*21*21 

- *12*12 

- ^22^22 
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In  order  for  this  condition  to  be  satisfied,  at  least  one 
of  the  following  equalities  must  hold: 


Otl^2  = Ct-2^1*  (21) 


For  m - 2 we  have  the  equations  (16),  the  equations 
obtained  for  m = 1,  and  the  following  additional  equations: 

^11M01  + ^21*n01  ~ ^12^02  ~ S22n02  - 6c2(S11^rt2i  + *21>N2l 
— - ^j_2^^22  ^22^^22 ) = 0 

^11^11^01  ♦ ^21621n01  - ^12612m02  ~ ^~22622N02  - 2°2 

*^115611M21  * ^215€21K21  “ <5125€12M22  “ 622^£22N22^  = 0 

£lilsilM01  * ^211-S21%1  - *121^12M02  - ^221^22N02  - 6°2 

* ^ ^11^11M21  + *21S21N21  * S12S12M22  “ ^22^22N22^  = 0< 

(22) 


Unless  (21)  is  satisfied,  we  must  tate  M11  “ N11  ~ M12 
= %2  = 0.  Substituting  (17)  in  (22)  gives 


- 18 


SllM01+  s21n01-612m02  = <522[N02+6c2^22-612)(cii1-a21)N2d 

^ll£llM01i‘<S21€21‘ n01-612€12  ^02  = 5 22^22^02 

+ 2c2(^22~€12)  ( ^ll2,f  ^212-  ^122"^222  + f5laI^’“^2ci2^  ) N22  ] 

5ll1SllM01+S2l1S21H01-^12ls 12M02  = S22 [S222S22N02 

+6c^(e- 22“  612^  ^lal^"/92a2^^22j  • (23) 


These  have  the  solution 


Mqi  — 


522^2-^1^1  + 
fillL€21-£ll 


(£21-£12HK  3^2-2^!) 
2(^2-^l)(£21-£ll)  - 


*01  = 


622 

K2"£12K1 

*>11 

€21-^11 

*02  = - 


2(V2->'i)  (£21-€ll) 

<522 

K 3 +K  2 - 2 ^ 

fi12 

2(*'2-yl>  J 

ft 

(24) 


where 

Kl  • Nq2  + 6c2(^22_612)  (ctl^“ct'2^^N22 

K.2  = €22N02 

♦ 2c2  (^22“£12^  ^112  + ^212~  ^122"^222‘f^la  \^~^2a'2‘  )N22 
K3  = ^222s22N02  + 6c2(<=22-£l2)(./3lal1-^2a21)N22. 


19  - 


Thus  we  see  that  the  only  polynomials  which  give 
exact  solutions  are 

f = PoU2-^2)  + QoU2-£22y2) 

F = - 2M0<5]_xy  - 2N0»2xy 


F = -AM2^1xy(x2-512y2)-2M0S1xy-4N2S2xy(x2-^22y2)-2N0^2xy* 

(25) 
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4.  approximate  solutions.  3ince  so  few  polyno- 
mials give  exact  solutions,  we  must  have  recourse  to  approx- 
imation. One  method  of  obtaining  approximate  8olutlona  la 
by  the  uae  of  function  space. ^ fie  will  give  a brief  intro- 
duction to  the  parte  of  this  method  which  apply  to  our 
problem. 

-4 

A vector  P in  the  function  space  is  the  ordered  set 
of  functions  (fxx,  Tyy,  -fXy)  which  represent  a state  of 
stress  of  the  body.  The  scalar  product  of  two  vectors  is 

given  by 

= */*(  e ii^n'  ♦ 2e^2^12'  + e22^22'  )dv 
~ ^®ll  '^ll  * 2®12 '^12  * e22  ' ^22  (1) 

where  the  ejj  are  connected  with  the  ^ijby  Koose's  law 
(1-5)  and  the  integration  is  over  the  entire  volume  of  the 
body.  The  length  of  the  vector  ? is  given  by 

P2  = |p|2  = P-P  (2) 

The  distance  between  two  vectors  Is  the  length  of  the 
vector  connecting  them.  Two  vectors  are  said  to  be  orthog 
anal  If  their  scalar  product  vanishes. 

^ V. . Prager  and  J.  L.  Synge,  ''approximation  in 
elasticity  based  on  the  concept  of  function  space." 

Suerterly  of  Applied  Mathematics  5.  241-69  (1947). 
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An  exact  solution  must  satisfy: 

(1)  the  equations  of  equilibrium, 

(2)  the  equation  of  compatibility, 

(3)  the  boundary  conditions,  and 

(4)  the  continuity  conditions  at  the  layers  of  separa- 

tion. 

The  stress  state  of  such  an  exact  solution  is  called  the 
natural  state  and  is  denoted  by  S. 

If  the  boundary  conditions  are  conditions  on  the 
stresses,  or  if  such  displacements  as  are  specified  are 
zero,  the  problem  is  said  to  be  one  with  stress  boundary 
conditions.  In  this  case  we  have  also  the  following  nota- 
tion: 

"§»  is  called  the  completely  associated  state.  It  sat- 
isfies the  equations  of  equilibrium,  the  boundary 
conditions  on  the  stresses,  and  the  condition  of 
continuity  on  the  stresses  at  the  layers  of  sepa- 
ration. 

Sp'  where  p = 1,  ...,  m are  called  the  homogeneous 
associated  states.  They  satisfy  the  equations 
of  equilibrium,  the  conditions  of  continuity  on 
the  stresses  at  the  layers  of  separation,  and 
have  zero  stresses  at  points  on  the  boundary  where 
the  stresses  are  specified. 
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S " where  q = 1,  . ..,  n are  called  complementary  states. 

They  satisfy  the  equation  of  compatibility,  the 
conditions  of  continuity  on  displacements  at  the 
layers  of  separation,  and  give  zero  displacements 
at  points  on  the  boundary  where  zero  displacements 

are  specified. 

tp'  where  p = 1 m are  called  orthonormal  homo- 

geneous associated  states.  They  are  obtained  from 
the  set  of  homogeneous  associated  states  by  ortho- 

normalizing.  -1- 

f " where  q = 1,  . . . , n are  called  orthonormal  comple- 
mentary states.  They  are  obtained  from  the  set  of 
complementary  states  by  orthonormalizing. 

It  is  easily  seen  that: 

I.  The  homogeneous  associated  8tates__f orm  a linear; 
ma  ni  fold  7H  called  the  homogeneous  associated  aanlfplo^ 

I i.  If  any  homogeneous  associated  manifold  b^_added 
to  a completely  asso^l  state,  the  result — 1 s ..a — Sjorq et.e  1 v 

r eclated  state.  -The  subapace.of  completely  associated 
Bt.at«a7n»  is  called  the  completely  agqoclateq.  sqbgPC.aS^- 

III.  The  complementary  states  form  a linear  manlfqlj_ 

called  the  complementary  manifold. 


1 A convenient  method  of  orthon'>r’®eJ^^G4n-6iVen 
in  U.  0.  Peach,  ^Simplified  technique  for  constructing 

orthonormal  functions."  Bui,  Bmer. Math . — ooc^_  ^0,  --'b-  9 
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The  following  theorems  are  proved  in  Frager-Synge: 

IV. The  homogeneo ub  as soclated  manifold  and  the 

complementary  manifold  are  mutually  orthogonal. 

V_. For  the  natural  state  ia  and  any  completely 

aasoclated  state  o». 


S-(S  - 3»)  = o.  (3) 
Hence  the  extremity  of  S Ilea  on  a sphere  having  S»  for  a 
diameter;  the  center  la  at  jro>  and  the  radius  la  ^8». 


VI 


The  natural  state  S la  orthogonal  to  any  homo- 


geneous  asaoclated  state  b 


o • o ' — 0, 


Lhl 


Hence  a la  orthogonal  to  the  homogeneous  as a oclated  manl- 
fold  7 7)  and 


3«Ip'  =0  (p  = 1,  m) . (5) 

VII.  The  difference  between  the  natural  state  S 
and  any  completely  aaaoQla.ta<3_atate .^»-lfl-artilaas2a6l-la 
any  complementary  state  ~S" 

(3  - 3»)-S"  = o.  if) 

Hence  (a  - 3»)  la  orthogonal  to  the  complementary  manifold 
71  and 

3-V1  = 3«-I<r  (q  = 1,  n).  (7) 

VIII.  The  difference  between  the  natural  atate  S and 
any  completely  associated  state  Ha*  la  orthogonal  to  the 
difference  between  a and  any  complementary  atate  o'1; 

(a  - ia*)*  (a  - ^" ) = 0,  (8) 


2b 


Thus  the  natural  state  lies  on  every  sphere  which  has  for 
Its  diameter  a line  joining  a point  of  the  associated  sub- 
space  Tfl * to  a point  of  the  complementary  manifold  7). 

Obviously,  the  uncertainty  Ln  the  position  of  3 
will  be  least  if  we  pictc  a point  of  Tfl*  and  of  71  such  that 
the  distance  between  them  is  a minimum.  To  do  that  we  let 


m 

M*  = 3*  ♦ E a I ' 
P-1  V 


be  a point  of  771*  and 


N 


n 

= E tqTq' 

q=l  ^ 


be  a point  of  71.  The  square  of  the  distance  between  them 
is 

(fi*  - N)2  = (3*  ♦ Z aplp'  - E Vq")2.  (9) 

P=1  q-1 

A necessary  condition  for  a minimum  is 

ap  = - 3*- Ip',  bq  = S*.Tq". 


Since  7m  and  7)  are  mutually  orthogonal,  this  is 
necessarily  a minimum.  V.'e  have  then 

v = i»  - r (s*-ip')ip' 


P-1 


n 


V = £ (s*.iq")V- 

q-i 


(10) 
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lie  define  the  error  a of  an  approximation 
3 - 5 . If  any  point  on  the  sphere  which  has  as 
eter ’ (Vm*  - vn")  be  chosen  os  3,  then  we  have 

«2-(v--V>2* 

If  weciiooss  for  3 the  center  of  the  sphere,  which 
bJ  0 = i(V.»  * Vn") 

the  error  satisfies  the  inequality 

»2  « i(v  - V)2  = i‘V2  ~ 7n"2>- 


3 by 

its  diam 

(ID 

is  glven 
(12) 


(13) 
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Application  to  laminated  construction.  A s a 

completely  associated  state  we  can  tase 


S*  = (V.  ryy*,  rxv») 


'xy' 


= (' 


V ax2* 


c^F»  \ 

dx^j  ) 


(1) 


where 

F*  = RfF^^)  + F2«(z2)} 

with 

Fl*  = Z (F n*  -►  i^n*^zin+2 
n=0  n i 

F2*  = Z (Qn*  + lNn^)z2n'2 

n=0 

and  the  Fn*,  Qn*,  ^n*>  an<3  satisfying  the  equations 

(3-1),  (3-2),  (3-6),  and  (3-7). 

% 

Similarly  for  our  first  complementary  state  we  can 

ta  tee 

3!"  = (W.  tjy" ' W)  = (*Sg.  2f^,  - &£)  <’> 

where  F-^"  is  defined  In  a manner  analogous  to  F*  except 

that  equations  (3-3),  (3-4),  (3-5),  (3-3),  (3-9),  end  (3-10) 
must  be  satisfied. 

Then 

4 

Vq*  - J* 


(3) 


Since 


?i"  = hm/  ft"i 

we  have 

V1 

S1  *S1 


(*) 


Using  equations  (1-19)  and  (1-20)  we  obtain 
S*«  ~ J1  ^lRl*Rl"  ■*  &]_2^ + 5^^ 

+ 82262R2*R2"  + ^lRl*Rl"  ♦ £2^1*^2"  + ^1  R2*Rl" 

+ ^2R2*R2"  + 2 V{  S-^I + £^S2^1*^2" 

+ ^1^2^2*^1  ' * ^2  ^2*^2  )J  (5) 


where 


Rl*  = the  real  part  of 


fr2' 

5z22 


<3^Fl* 

aZl2 


R2  = the  real  part  of 


l2"  = the  imaginary  part  of 


^2  7 .1 

5 J2 
bz22 


Using  (1-22),  (1-24),  and  (1-11)  equation  (5)  can  be  written 

= ^[^(R^Ry'  + Ii»tirt)  + 8'22(R2^R2"  ♦ 12**2") 

* (i  + i^)(Rl#R2"  + R2*Rl")  ♦ 5!52(*1#*2M  + *2*ll")[dv* 

(6) 


In  a similar  manner  we  obtain 
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3*  * - 

j2v|Sy^(Ri#^  + ty*2)  + + t2*^) 

♦ 2(i  * f£K**2*  * 2SlS2t1’I2*]<iv 

(7) 

and 

< l<  *» 

31  ®1 

= JZV[S l2(Hl"2  ♦ il"2)  - ^22(h2"2  * J2"2) 

* 2{£  *■  yxjl<l"R2"  * 2S1S2l1"S2"J  <3v 

(8) 

If  we  tatce 

s - HV0*  ♦ V!"), 


the  maximum  error  can  be  found  from 


«2  ^ k 


S* • 3*  - 


(1»-V)2' 

. 3 " * 

3i  3i 


(9) 


i»ny  set  of  Fp'  such  that  Tpxx' 


a^Fp' 

a v2 


rpyy 


' = 172“'  snd  T> 


pxy 


I _ 


>2  -p  i 

o *p 
^xay 


equal  zero  on  the 


boundary  and  such  that  fpyy ' and  'fpxy'  are  continuous  at 
y = ic  will  serve  to  determine  the  *p ' . Such  a set  is 
given  by 


Fb’  = (x2  - a2)2(y2  - b2)2f „( x , y ) 


(10) 


where  fp(x,y)  Is  continuous  and  hae  continuous  derivatives 
in  the  closed  region  occupied  by  the  body. 

i*  set  of  io  ' may  be  constructed  from  any  set  of 
displacements  which  tane  on  zero  values  at  points  of  the 
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boundary  where  the  displacements  are  specified  to  be  zero, 
and  which  are  continuous  at  y = *_c. 


It  is  not  necessary  that  the  F*  used  to  determine 
the  completely  associated  stated*  satisfy  equation  (1-15). 

ryy 


2 

any  F*  such  that  'f’xx*  = d 


'*  = and  Txy* 


a2F* 


dy^ k JJ  2>x2>y 

satisfy  the  boundary  conditions  and  such  that  Tyy*  and  TXy* 


are  continuous  at  y = +c  can  be  used  to  obtain  an  approxi- 


mate solution. 


(L Illustrative  examples,  Let 


*1  = ■i?lzl'5  f2  = ^lz2"  • 


(1) 


Then 

3*  = [-6(S]2Px*-*-^^i» )xf  6(Pi*  6( Si2Pi«-+  ^22Qi* )y]  (2) 

where  the  values  of  P^  and  In  the  first  layer  are  related 
to  their  values  in  the  second  layer  by  the  equations 


*11*  * all”  - *12*  * «12*  = 0 
11 # * ^21^^11*  ” ^12^^12*  “ ^22^12*  = 0.  (3) 

These  have  the  solutions 

pH*  = [(  ^212- ^122)-t’l2*  + (&212"  i222)^12^]/^212-<5ll2) 

= [(S^2^~^11^)F12#  + (^22^“  ^11^  )^12^^21^_  ^11^  ) * 

Also 


- [-6(  S12P^"-*-S22Q]_''  )x,  6(P^"-*-Q^  ' )x,  6(S;l^Pi  +^22Qi')y] 

(5) 

with  the  equations 

ellpll"  + - €12f12"  - e22^12''  = 0 

SllFll"  + s21^1l"  - S12P12"  - 522^12"  = °»  (6) 


which  have  the  solutions 
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*ll" 


j.€;2S21"g215 12)^12"  + (C22S21-*21S22)QI2  ' 


^ 5 — 6 5 

ir2i  21  11 


'*‘11  - 


^€llS12-el25ll)P12’'  + ^*11S22~  €22S11)Q12  ' . 


(7) 


C11S21 


e21Sll 


Let  ua  further  require  that 

P-^2*  Q]_2*  ~ ^ 

6^2^i>i2*  * *22^^12*  = ® 


and  ta£e 

Pl2"  = P12*  and  Q12"  = *12* 


(8) 

(9) 


These  £ive 


p12 

" = *12* 

' = ^222k-/(  ^222- *122) 

^12 

" = q12* 

► = -6122K/(6 

222-'5122 

) 

Pll*  = 

^212V  (^212- 

Al2) 

Qll*  - - 

• &112K/  ( $21~- 

(1C) 

txxl#  = 

0 

ryyi* 

= 6Kx 

Txyl*  = 

0 

exxl*  = 

-6pi Kx 

eyyi* 

- 67'iKx 

exyl*  “ 

0 

Txx2*  = 

0 

Tyy2* 

= 6xvx 

7^2*  = 

0 

exx2*  = 

-6/S2^x 

eyy2* 

= 6Y2&X 

e xy2*  = 

0 

(11) 

- 32  - 


«r  " - 
'lxxl  ~ 

6(^2-^^  )Kx 

^1^1  “^1 2 

elxxl '*  - -6/^Kx 

1'lxx2"  = 0 
elxx2  = -6^2^-x 

= ^]qX2  [io56,r2^2dy+ Jp36r2^2dyjdx  = 48Y'2K2a>t  (13) 

3**S*  = ^ )qX2  [j^36triK2dyr</^36r’2K2dyjdx  = 48ri2a^c(  Y\-  *£)■♦■  l^b] 

(1*) 

V'V  = •,*Io(lo[361i2x2t'S22^-2/3l/S2^^2':(l)/(al^l-%2) 

f 72  t^li2y2(|S2>r’l-'5]/2)2/(a]/ i-^2)2]iy>Jj’6r2K2x2dy)jjt 
= 48K2{  [(/«22r1-^iy327'2-a1/22)/(a1/1-/S12)]a3c 
♦ [2J/1(^1y'2-^2«'1)2/«»l'21-d12)2]ac3  ♦ Y|o3(b-o)}.(15) 

As  a numerical  example  let  us  consider  spruce  with 
the  longitudinal  axis  in  the  direction  of  the  x-axis  in  the 
middle  layer  and  the  transverse  axis  in  the  direction  of 
the  x-axis  in  the  outer  layers,  the  radial  axis  being  in 
the  direction  of  the  y-axis  in  both  layers.  Let  us  suppose 
that  the  layers  are  of  equal  thickness  and  that  the  width 
of  the  piece  is  ten  times  its  thickness.  Tie  assume  the 


'lyyi 


•lxyl 


»«  _ 


Qiri'^L2  al^l‘^l2 


elyyl  ~ 6v^Kx  elxyl  ~ Vl  ^lxyl  ' 


^yy2  ~ ^Kx 


^lxy2 


,q"  = 0 


el  yy  2 " = 6f2^x  elxy2".=  0 


(12) 
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following 

to  be  the 

constants  of 

the  material 

M 

rH 

0 

0.507 

•10"6 

d2  = 

14.3-10"6 

Z3!  = 

0.314 

•10'6 

@2  ~ 

4.48-10”6 

*1  = 

6.48* 

lO-6 

y2  = 

7.69-10*6 

II 

4.81- 

10-6 

v2  = 

108.7-10-6 

Vie  calculate  the  following  results: 

= 29.9c4K2,  2* -3*  = 28.3cAft2!  51”-31"  = 78.6cAK2 

Vx"  = 0.38o3x",  « - 2.05c2ii,  a/area  ^ 0.006K. 

We  obtain,  for  our  first  approximation: 

In  the  first  layer 

Txxl  = -9.50Kx  -fyyx  = 5.35KX  <xyl  = 9-50Ky. 

In  the  second  layer 

rxx2  = 0 Tyy2  = rXy2  = °- 

If  we  wish  an  approximation  which  satisfies  the 
boundary  values  we  taxe  as  our  solution.  This  gives, 
for  both  layers 

'T'xx  = 0 ^yy  = 6Kx  ^xy  = 0 

with  an  error  of 

e/area  - O.OIK. 

To  improve  this  approximation  we  tafce 
Fx’  = (x2-a2)2(y2-b2)2x 


This  gives 
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■§1’  = [4x(x2-a2)2(3y2-b2) , 4x(5x2-3a2)(y2-b2)2, 

-4y(x2-a2) (5x2-a2)(y2-b2)] 

= 256  [(32/3465)c‘3a11(9c5/5  - 2b2c3  * b4c) 

+ ( 16/315) ^a9(3cV7  - l^cb/5  + 5bV/3  - b6c) 

*■  ( 1/7)  V5 a?  ( c9/9  - 4b2c7/7  + 6b4c5/5  - 4b6c5/3  ♦ b8c) 

+ (64/315)^5a9(c7/7  - 2b2C5/5  + bV/3) 

+ (123/17325)02^5  _ ( 512/33075 )/32a9b7 
+ (123/2205)  r2a7b9  + (5l2/33075)^2a9b?] 

where 

a3  = ax-d2,  ^ =^1-/S2,  ^ = )<1- ^2,  = • 

3*  • 3X  ' = -(  256/35 )/®3 a7 ( c^-b2c)K 

1 1 ’ = sl'  / 1^1  'I 

= s»  - CS*-t1,)i1'  = - (3*-31,/l1' •31' )^x' . 

For  the  numerical  ccae 

31',3,1'  = 9.93-1015c16  3**^'  = -5.34*106c10K 

3i'  = [l2x(x2-900c2)2(y2-3c2) , 20x( x2-540c2 ) ( y2-9c2 ) 2 , 

-20y(x2-900c2)(x2-l80c2)  (y2-9c2)] 

V = [4.69(x/30c)(x2/900c2  - l)2(y2/3c2  - l)cK, 

180  + 1.41(x2/540c2  - l)(y2/9c2  - l)2  (x/30c)cK, 
-4.70(y/3c)  (x2/900c2  - 1)(x2/130c2  - l)(y2/9c2  - 1)ck] 
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As  a second  example,  let  us  consider  the  same  prob- 
lem except  that  the  boundary  conditions  are 

for  x = +30c  (30c=a)  1^x  = K(1  - y2/t2),  <Xy  = 0 

for  y = +3c  (3c=b)  Tyy  = 0 rxy  = °* 

These  boundary  conditions  will  be  satisfied  if  we 

taKe  F*  = |-K.y2(l  - y2/6b2). 


If  we  also  taice 


(x2-a2)2(y2-b2)2 


F2'  = x2(x2-a2)2(y2-b2)2 
F, ' = y2(x2-a2)2(y2-b2)2 


we  obtain 


3*  = [K(  1 - y2/b2) , 0,  o] 


J1*  = [4(x2-a2)2(3y2-b2),  A^-a2)^2-!*2)2, 

-l6xy(x2-a2)  (y2-b2)] 


S2'  = |>x2(x2-a2)2(3y2-b2) , 2(15x4-12a2x2+a4)(y2-b2)2, 

-8xy(3x4-4a2x2+a4)(y2-b2)] 

53'  = [2(x2-a2)2(15y4-12b2y2+b4),  4y2(3x2-a2)  (y2-b2)2, 


-8xy (x2-a2) (3y^-b2) (y‘ 

= -2.86*  104c8K. 

3i' *§i’  = 1.02-1012c14 

S*.^2*  = -3 . 68 • lO^cb^A 

tx' •^2’  = 7.89-10l3cl6 

= -6.59*104c10K 

‘31’  = -5.90-1010c16 
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22'-V  = 1.68*10l6c18  %'•%'  = A.38-101M8 

= -7.06*10l2cl8 

Orthonormallzing  gives 

I1'  = 9. 90*10-7c-7?  ’ 

I2'  = -7.48*  lO-^c-T^'  ♦ 9.67*10-9c-9^2' 

I5'  = 1.90*  10“9C"7^1  • + 1.12-10-11c-932'  ♦ 4. 77*  10-8c*93j ' 
Then 

= 3*  ♦ 1.74*10"8c-6K^1'  + 1.37*  10_10c“8fc.ii2 ' 

+ 1.55*  10“10c“8K^  ' . 

We  also  tase 

Si"  = (o>,  0) 

then 

= 2.58-1C)8c2  = 240Kc2 

tj/'  = 6.23-10-70-1^"  Vj/'  = 9.34*10-7kI1" 

rixxl"  = I-**.  Tlyyl"  = 0.092K,  rlxylM  = 0 

'rixx2"  = 0.05534,  rlyy2"  = 0.03224,  'flxy2"  = 0 

4s  a third  example  let  us  consider  the  same  piece 
of  spruce  with  the  boundary  conditions 


for 

x = 0 

Txx  ~ o, 

Txy 

= A(1 

x = 60c 

u = 0, 

v = 

0 

for 

y = T.  b ( b=3c) 

T =0 

^yy 

^xy 

= 0 

These  are  satisfied  If  vie  taite 


F*  = -Axy(l  - y2/3b2) 


giving 

S*  = [2Axy/b2,  0,  A(  1 - y2/b2)]. 
Tie  further  tate 


Fx'  = x2y(y2  - b2)2 
?2 ' = x^yly2-  b2)2 
' = x2y^(y2  - b2)2. 


These 

give 

t • 
°1 

[4x2y(5y2-3b2),  2y(y2- 

t2)2 

, -2x(5y2-b2 

) ( y 2 - b‘ 

h' 

= 

[4x5y( 5y2-3b2) , 6xy(y2 

-t2) 

2 _• 

3x2(5y2- 

b2) (y2 

% 1 
s3 

[2x2y( 21y4-20b2y2+3b4) 

. 2y 

5<y2 

-b2)  , 

-2xy2(7y2-3b2 

• 

-> 

Si 

’ = 7.30  -102c7A 

? i 

bl 

•Sl’ 

= 7.62* 

107c12 

-♦ 

3*  • 

s2 

' = 3.43‘10ac8A 

3 1 
51 

•12’ 

= 3.80* 

10* 

-¥ 

s** 

h 

• =—4.28*  lO^c^A 

ai 

•V 

= 3.  6l- 

108c14 

? 1 

a2 

•3. 

2'  = 1,95  •1011c14 

q 1 

s3 

•V 

= 2.71  • 

109c16 

2 • 
°2 

•I- 

5’  = 1.83*1010c15 
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?x  ' = 1 . 15  • 10“'4c_631  ' 

12'  = -6.73  • 10"  V6^’  + 1.35  •lO-^c-^' 
ly  - -7.31  •l0~^c“^1'  - 2.92 • IO‘^c'^32'  + 3.22 -1Ct5c- -'£>2* 

V3*  = 3*  - 3-02 -10-50-5^’  - 3.32-10“7c-6A32' 

♦ 7 -82  • 10“^  ' 

The  a j " must  satisfy  the  compatibility  equation 
(1-4),  give  continuous  displacements  at  J =±c,  and  satisiy 
u=0,  v=0  at  x=60c.  One  way  to  guarantee  this  is  to  start 
with  a set  of  suitable  displacements.  Such  a set  is 

u = (60c-x)f(x)g(y) 
v = ( 60c-x)m(x)n(y) 


where  f(x),  m(x),  g(y),  n(y)  are  bounded  and  continuous. 
Let 

Ul"  = oOc-x  y-i”  - 60c-x 

u2"  = (60c-x)y  v2"  = ( 60c-x)y 


These  give 

ft  _ i 

elxx  “ 

eiyy 

^lxx  ~ 

riyy 

e2xx n = -y 

e2yy 

r2xx"  = -“>y+ 9*(60c-x) , 

r2yy 

0 0lxy"  ~ 

* / Tlxy  = _/*xy 

( 60c-x ) e2xy"  = |-(  60c-x-y ) 
-^y  + X(60c-x) , 

= V600-*-*) 
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§l"  - {-CO,  -j,  -/fry) 

a 2 " = [-u>y^(60c-x) , -^y-»-A(60c-x) , /<xy ( 60c-x-y ) J 


2*.V  = -2.40*102^C2  V*^1M  = 2.71*  10®c2 

3**32"  = -1.44*104Ac3  o-l"*S2"  = -2.77*109c5 

S2" -32"  = 7.12-I0l°c4 

‘l1M  = 6.08-10-4c-1S1'' 

I2”  = 4.95‘10-5c-1S1"  ♦ 4.32-10-6c-232" 

^2"  = -9. 64-10-5^^"  - 3.35*10-7Ac-132''. 

For  the  solution  we  taKe  the  stress  state  given  by 


C = i(V^  + V2"). 


The  values  of  the  parameters 

In  the  middle  la^er 

«>=  2.03*106 
<f>  = 9.85-104 
A = 1.59-105 
/Uxy  = 1 . 04*  105 


In  Si"  and  o2"  are 

In  the  outer  layers 

• «;  = 5-92-104 
<(>  - 3 • 45*  104 
A = 1 . 10* io5 
Mxy  - 4. 60* lo5 
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7.  Approximation  by  finite  equations.  Another 
method  of  approximation  Is  to  replace  the  dl 1 f erentlal 
equation  by  a set  of  finite  equations.  To  do  this  we  first 
approximate  F(x,y)  by  a polynomial  G(x,y)  such  that 


GUi.yj)  = f(xi»yj)  I1) 

at  certain  points  ( x ^ , y j ) . 

The  given  points  are  numbered  in  such  a way  that 
if  the  lines  x = x^  and  y = yj  for  1 = 1,  I end 

j = i,  J are  drawn  In  the  xy -plane,  the  point  (xi,yj) 

is  the  intersection  of  the  lines  x = x^  and  y = yj.  Then, 
from  Lagrange's  formula  of  interpolation,  the  required 
polynomial  is 


GU.J)  =zzu 


n'  (y-yS2)  n'(x-xi»2) 


r2 


8lrl 


n'(y8l-ya2)  n'  (xri-xr2] 


T(xr,ya)  (2) 


s2 


r2 


where  IT'  indicates  that  in  the  product,  r,  is  not  to  to*e 

rJ 

any  of*  the  values  rj  = r^,  r2 rj_i.  If  ® is  the  only 

value  in  the  range  of  r2,  we  tatfe 

rr ' (x-xr2) 


n ( x — x ) 

Similar  expressions  Involving  y and  s have  similar  meanings. 
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That  G(xlfyj)  = F(x1(yj)  ccn  be  3een  froin  the  fact 
that,  If  x = xm  and  y = yn,  where  m belongs  to  the  range  of 
r and  1 ^ n ^ J,  then  the  following  hold: 


n'(x-xr  ) 

r2  _ (0  fsr  rir  ® 

n'(xri-xr2)  (1  for  rx=  m 


n'  (y-y8  ) 

82 

n (ygi-ya2) 


(0  for  ai^n 
(1  for  8}_=n 


V. e replace  the  derivatives  of  F in  the  differential 
equations  by  the  derivatives  of  G given  by  the  following: 


n'ly-y  ) rn'U-x  ) 

30  = Z Z ^ — ^ L F(xr  ya  > 

H 3iri  n’(jBl-ya2)n'(xri-xr2) 


(4) 


32G 

dx2 


= 11 


n'(y-yB  Jl'Il’ri'U-x  ) 

a _ tv*  r-  r..  ^ 


r2  r3  r* 


!r!  n'(ySi-yS2)n'(xri-xr 


F ( X > y 8^^  (5) 


r2 


«^G 

3y 


5*=  II 


r,n'(y-yS5)ri,(x-xr2 


) 


b2  8; 


a1r1n ' (ySl-y82)n ' (xri-xr2) 


F ( x , y8l) 


(6) 


42 


32G 


= HZ 


X ' X ' n ' ( y-y3 , „ 

s } 4 r 


)n ' (x-xr2: 


8lrl 


* ( *r-^ » y ^ 


(7) 


_a_^G  a 

dxdy 


IZ4I3 


r'n'(y-yB  )Z  'n’(x-xr 

s ^.3-2  u r~  r,  -> 


) 


2 *3 


3lrl 


n'(yBl-J32>n'(xri-xr2) 


xr^»  y 8-^  ^ 


(8) 


Let  the  lines  x=x^  end  y=yj  be  evenly  spaced  at  a 
distance  of  d apart,  and  let  the  values  of  F^j  = F(x^,yj) 
be  given  at  the  points  shown  below. 


'-1.1 

0,2 

Jo.i 

Fi,i 

F1.0  . g 

o 

CM 

1 

1*4 

F-1.0 

*0,0  _ 

F0^1_ 

Fl,-1 

F 

*0,-2 

(9) 
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Using  (5),  (7),  and  (8)  we  obtain 

=(1/I2d2)  [-(exx)0)_2+l6(exx)0t_1-30(exx)0>0 
d y^  / o,  0 -« 

t-l6(  exx  ) o , 1"^  exx ) 0,  2J 


(10) 


^®-p)  =(1/I2c2)  [_-(eyy).2,o+l6(eyy)-l,0->0^eyy)0,0 


V a*2 


0,0 


+ l6(  eyy )l , 0"^  eyy ^ 2,  o] 


(ID 


(Sfio  o=u/3d2) 

’ (12) 


Then  using  (1-13) » (5).  (7),  and  (8),  and  substi- 
tuting the  results  in  (1-4),  we  get  the  equation 

<-l44^0>_2+96a0>.1+900a0>  0+96^^-144^^-1300,00^-1447:2,0 
'■96'C1,01-900r'0jo'-96r'ij0-14Ar8j0»216(y.li.1*^1|.i»^iii^ltl)>Fo>o 
* (480^lj0,*80^0|0^10AY-,2_0-320r1>0-43or0i0t64r'lj0>56r-2>0 
-120V_1>_1-72V1_1-120^_1>1-72^j1)F.1>0  + ( 1040^  _2-320c<0,  _x 
-48000^  0+ 64^ ^ 1+56a0> 2+384^0>  _i-48^_2,  o*480/^,  0"48/S2, 0 

-160v_1>.i-160^>.1-96z:i>i-96V1>i)F0>.i  + ( 480/?0>  0+480£l,  0 

^•56r_2, 0+64Yli,  0-480^0, 0-320ri,o+104r2, 0-72 Z/.1,_1-120^1,-1 
-72y-l,l-l20l/l,l)Fl,0  + (56a0)_2.64aQj_1-430a0>0-320a0>1 
+l°4c'0,2f584^0>1-43/8_2>0+480/90>0-43^2>0-96V_1>1-96z71>_1 

-16'OV-i,  1-160^1, i)F0,i  ♦ (-192^0>. i+4^_ 2,  o-256/3.i,  o+1^2fO 
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+ 120Z''_1  ( 40n,-l  + '72^_1,1+24^>i)F.1,-i  ♦ (-192/0,-1 
+ 16/_2,0-236/i,0^^/62,0  + ^-1,-1  + 120i/1.-1  + 24^*1+72^1':l)F1’'1 
i-  (-192^>1^l6/B_2,o-256/l,0+W2,01-24V.i,_i^72^i)_i  + 40n,l 
*120*'l,l>?l,l  + (-19^0,l  + 48^-2,0-256/-l,0+l6/A2,0+^-l,-l 

+24^  _1i'1209_1>1-»-4091>1)F_1>1  «■  (-30/0>0-30/4.2,0-’5/-2,0 
♦176^,0^0/b,  0-16^1, 0-1^2, 0-36v.i,.1n2^1>.1-36^1,i 

+ 12^j1)F_2>0  -+  (-35a0,-2i‘1’76flt0,-l^0oi0,0-l6cl0,l-llc‘0,2 

-/S-2,0+l6/S-l,0-50A0,0  + l6/fil,0"^2,0)F0,-2  + ( -30/?0,  0-30/?2, 0 

-11/.2, 0-l6^-l,0r50/0,0  + l76/l,  0-33^2, 0 + 1^-l,-l-56l/l.-l 
♦12t/_i,i-36Vi,l)F2,0  * (-Hao,-2-l6Q(0,-l+30^o,0+136«0,l 

-35o»o,  2“A-2,  O^-l^A-l,  0-30/0,  0+l6Al,  0-A>,  2)fo,  2 = °-  (13) 

If  a,  p,  i',  fcnd  v are  Independent  of  the  first 
subscripts , we  have,  where  the  second  subscript  Is  that  given 

(-I44a_2*96a_i  ^900ao+96a1-l44a2-i800/s0+304/o+432V-1  + 4323>L)Fo,0 
+ (960/o-576/q-192»/.i-192i^)F_1>o  ♦ (104a.2-320a.1-4Boc»0 

>64a1i.56tf2^V-1+334A0-320^1-192>/1)F0>.1-*-  <960/90-57 6r0 
-192V_1-192^l)F1>0  + (56a<.2»640».1-480«0-320«1+104‘»2O84^0 
334/51-192i/.1-320^)F0>i  - (-38^30+1602/_i^96^i)F-i,.i 

* ( -364^0^160  U.^96  )Flf.!  ♦ ( -384/J0+96z/_i+l60i'L ) F1 , l 

* (-364/b0«-96v_1+l60V1)F_1>1  + (-60/50>l44)r0-24J/_1-24i/1)F_2>0 

* ( -33a_2  fl76a_i+’30ot0-l6ci1-ll<i2 ) ^0,  -2  + ( -60/0^l44K'0-24z7_1 
-24Vl)F2,0  + (-Ha.2-l6o'-1^0«o>l76Dti-35ci2)Fo,2  = ^lA) 
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If  ell  of  the  points  lie  in  the  some  region,  we  have 

(804a-l800/a+804f+864i/)F0f0  + ( 96Q./S-576r'-3-84'T^r_1 , 0 
t-  ( -57&a<-768^-5l2v')Fo, -1  + ( 960£-576r-334vO Fi, o + (-576a 
♦ 763(3-512^)  F0>  x + (-364/?+256i/)F-1>  _i  + ( -334£+256i>Fi  , _i 
r (-584^+256^)?!^  + (-384/»i-256v)F.lfl  + ( -60/i+l44 /- 48  j/) F_2,0 
■*■  I44qr0  _2  *■  ( -60^+144 /■-48v')r2>  o + 144clFo,2  = 0 (*5) 

For  points  which  lie  on  the  boundary  of  two  regions 
we  ts <te  the  average  values  for  the  values  oi  a,  /3,  V,  and  V, 

If  10,  2)  lies  on  the  boundary  and  the  rest  of  the 

points  lie  in  region  I,  we  have 

(-72a[r  ■*-876aI-l300^x +304^x^3641^  )F0,o  ♦ ( 960/?r-576/i -334i^  )F-1,0 
+ (28arI-o04dI+768/SI+512^[  )F0,-i  ♦ ( 960ft -57o/i -334Vi  )Fi,  o 
, (52aII-o2ao<I+768^r-5l2VI)ro,i  + (-334^1 +256Vi)F_i,_i 
+ ( -384ft+25oft  )Fi, -l  ♦ (-334ft+256Vr  )b'i,  l *•  ( -334ft +256ft  )*  -1 , 1 
+ ( -60ft  +144/i -48vj  )F_2,0  + ( -5>Vft 1 +149^1 ) F 0,  -2  + (*^I 
1 1 44Yx - 48  ft ) F 2 , 0 ♦ ( -17^X1  +l6l£ft  )f0, 2 = 0 (13) 

If  the  point  (0,  2)  lies  in  the  region  II,  the 
points  (-1,  1),  (0,  1),  and  (1,  1)  on  the  boundary,  and 
the  rest  of  the  points  in  region  I,  we  have 
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(>96q(ii+900Q(i -l800/?£+804yx+2l6^xx  + 648z/I  )F0,0  + (960ySi -576% 
-96^II-238>/I  )F_1>0  + (88otII-664aI+T68/3I-96^II-4l6^I)F0>_1 

♦ ( 9^0^-576^ -96  VII-238VI)Flj0  + ( -56a11  -520*% +192^  j 
+ 576/% -160^i  £-352%)F0>1  + (-384%;+48%i+203%[)F_if  _i 

+ (-384/3I-r48vII+208vI)F1>_1  ♦ ( -384^+80*^ +176^  JF^ 

+ ( -334^+80^  r + 1761%  + (-60/eI+l44/I-12VII-36yi  )F-2>0 

( -190^+1630%  )?0>-2  * ( -60^+144)^.-121^-36^  )F2(  0 

♦ (53<»II*91<*I)F0,2  = °-  (17) 

If  the  points  (0,  2),  (-1,  1),  (0»  1)»  (1»  1 ) 

lie  in  region  II,  the  points  (-2,  0),  (-1,  0),  (0,  0), 

(1,  0),  end  (2,  0)  on  the  boundary,  end  the  rest  In  region  I, 

we  have 

( 403%  j +402%  -900%  j -900%  ♦ 402%  j + 402%  + 432^+432%  )FQ>0 
+ (430/^i  i-480/?x -238 Yu  +288/i-192j%i-192j/i  )P.lfQ  + (-120aXI 
-4560ix+192^h  + 576^i-192^x  x-320  Vx)Fo,-i  + ( 480%££ +430/% 
-288y'iI-233/x-192iAI-192^t)Flf0  + ( -456% : -120%+192/% j 
+ 376i%-320^[i-192Vx  )?0>1  + (-192/%i-192/%+92t%i+160%  JF^  ^ 
t ( -192^x1 -192/% +96^x +160^  + ( -192/%  r -192% +160%I 

+92%)F1>;L  + ( -192^x -192^+160%  i+96%)F_1>1  + (-30/?n 

-30/sI+72fII  +72%-24%£-24%  )f-2,0  + ( -I2ax x +156%  ) FQ|  -2 

+ (-30^ix-30/?x^72/ii+72%-24%I-24%)F2f0  + (156^! 

-l2ax)F0>2  - 0.  (l3> 

The  other  equations  can  be  obtained  by  interchanging 


subscripts . 
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The  above  equations  can  be  conveniently  represented 
by  diagrams  similar  to  (9).  The  equations  are  then  obtained 
by  multiplying  the  value  of  F at  each  point  by  the  number 
Just  to  the  right  of  the  point  on  the  diagram,  adding  the 
results,  and  equating  the  sum  to  zero.  The  following  Is 
an  example  of  such  a diagram  for  equation  (17)  using  the 
constants  given  on  page  33. 


027 


Region  2 

311 

-650 

311 

— ■ 

Region  1 

019 

-503 

1000 

_ * 

-505 

— 

1 

1201 

,-272 

201  - 

! -006 

We  calculate  the  value  of  F and  Its  first  derlva- 
tlves  by  Integrating  the  values  of  the  stresses  along  the 
boundary.  This  gives  us  our  boundary  values  for  F.  The 
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region  Is  then  divided  Into  squares  of  length  d by  the 
lines  x-x^  and  y-yj.  The  intersections  of  these  lines  are 
called  nodes . Nodes  falling  on  the  boundary  ere  called 
boundary  nodes.  If  some  of  the  lines  intersect  the  boundary 
in  points  other  than  nodes  we  cell  the  node  lying  on  the 
line,  interior  to  the  boundary,  and  closest  to  it,  a boundary 
node . The  value  of  F at  such  a boundary  node  la  approximated 
from  the  values  of  F and  its  first  and  second  derivatives 
at  the  boundary.  The  nodes  adjacent  to  the  boundary  nodes 
and  exterior  to  the  boundary  are  called  adjacent  nodes. 

The  values  of  F at  adjacent  nodes  ere  approximated  from  the 
values  of  F and  its  first  and  second  derivatives  at  the 
boundary.  Modes  interior  to  the  boundary  which  are  not 
boundary  nodes  are  called  interior  nodes,  an  equation  can 
then  be  obtained  for  each  interior  node  by  use  of  the 
appropriate  equation  from  (15)  - (18).  Thus  we  have  re- 
placed the  differential  equation  by  the  set  of  finite  equa- 
tions : 

1^jaiJF(Xl  ) = °‘  U9) 

In  general  the  number  of  equations  given  by  (19) 
is  too  large  to  maite  a direct  solution  practical.  The 
solution  can,  however,  be  approximated  by  the  so-called 


relaxation  method. 
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If  we  estimate  the  value  of  F as  F°  (taking  the 
value  of  F as  F°  for  a boundary  node  or  an  adjacent  node) 
and  substitute  F°  for  F In  the  left  sice  of  (19)  we  get 

Z »f^°Ui,yj)  = Ral  (20) 

If  R*1  = 0 for  all  c,  1,  then  F°  = F and  the  estimated 
solution  Is  the  correct  solution.  If  this  is  not  the  case 
let 


F’(xK.,y1)  = F°(xtt,y1) 


t Kl/n  Cl 


R'  v a 


HI* 


(21) 


The  cl’th  equation  then  becomes 


a^'Uic.n)  ♦ $K6ir0(xl,yJ) 


= R*1  - K'  R"L  - R 


Kl  _ u"Cl 


The  rest  of  the  equations  which  involve  F(xtc,yi)  become 


a^FUxi.yx)  - ^Kai  JF°(xi’ * j) 


rs  rs.-.icl/  cl 
R - aKlh  /&ZI 


Thus  by  changing  the  estimated  value  of  F(xt,  y^)  from 
fO(Xt>yi)  to  f 'lx^.y^  the ’Vesldual'  RKl  has  been  reduced 
to  h” ^ Starting  with  the  largest,  jhl,  each  is  relaxed  In 
this  manner  until  the  largest  remaining  ihl  Is  negligible. 
The  resulting  F’s  give  the  approximate  solution. 
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